Influence of the vacuum interface on the charge distribution in V 2 Oa thin films 
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The electronic structure of V2O3 thin films is studied by means of the augmented spherical wave 
method as based on density functional theory and the local density approximation. We establish 
that effects of charge redistribution, induced by the vacuum interface, in such films are restricted 
to a very narrow surface layer of « 15 A thickness. As a consequence, charge redistribution can be 
ruled out as a source of the extraordinary thickness-dependence of the metal-insulator transition 
observed in V2O3 thin films of ~ 100-1000 A thickness. 
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I. INTRODUCTION 

Over the last decades, V2O3 has been investigated ex- 
tensively by both experiment and theory - and still is 
a topic of great interest (see, e.g., Refs. [H, 0] and the 
literature given therein). Today, V2O3 thin films are re- 
ceiving growing attention, particularly because the thin 
film geometry appears to strongly modify the electronic 
properties of the material as is most clearly demonstrated 
by the famous metal-insulator transition observed in bulk 
V9O3 as a function of temperature, pressure or doping 
[3| . It is known that even a small alteration of the crystal 
structure, by Cr substitution Q or application of exter- 
nal pressure can have a serious effect on the metal- 
insulator transition in the bulk material. Therefore, it 
would not at all be surprising if structural modifications 
induced by the thin film geometry had similar severe im- 
plications. 

Recently, the observation of a thickness-dependent 
metal-insulator transition in ultrathin V2O3 films was 
attributed to the increasing c/a lattice parameter ratio 
(decreasing a, increasing c) because of interaction with 
the substrate [f|. However, this conjecture has been re- 
futed by both x-ray diffraction studies for high quality 
films @, IU and density functional calculations @- In 
fact, the experiments point to a very complex thickness- 
dependence of the lattice parameters with unsystematic 
changes of the c lattice constant for film thicknesses of ~ 
100-1000 A, thereby calling for a more in-depth study of 
the electronic properties of these films. Indeed, in other 
systems, like the LaAlOs/SrTiOs heterostructure [l(|, 
it has been shown that interfaces can seriously modify 
the electron density profile. A strong dependence of the 
electronic properties on the layer thickness has been also 
reported for the hexagonal layer compounds such as WS2 
jll| . In the present case a systematic investigation is still 
missing. 

In this paper, we aim at clarifying the influence of the 
vacuum interface on the electronic states of V2O3 thin 
films. At this interface, two possible scenarii may cause 



a significant reorganization of the electronic states: (1) 
a local modification of the crystal structure and (2) a 
charge redistribution due to the broken symmetry at the 
vacuum interface. In general, charge redistribution due 
to a surface, similar to heterointerfaces, may induce a 
considerable deviation from the bulk properties, see, e.g., 
Refs. EIOJ]. Even thou g h possibility (2) there- 

fore seems to be promising, it is the aim of this paper 
to point at the limitations of such a mechanism for film 
thicknesses exceeding some 15 A. 



II. COMPUTATIONAL METHOD AND 
TECHNICAL DETAILS 

The calculations are based on density-functional the- 
ory and the local density approximation (LDA). They 
were performed using the scalar-relativistic implementa- 
tion of the standard augmented spherical wave (ASW) 
method (see Ref. [l|| Chap. 2 of Ref. [l?] and references 
therein). In the ASW method, the wave function is 
expanded in atom-centered augmented spherical waves, 
which arc Hankel functions and numerical solutions of 
Schrodinger's equation, respectively, outside and inside 
the so-called augmentation spheres. In order to optimize 
the basis set, additional augmented spherical waves were 
placed at carefully selected interstitial sites. The choice 
of these sites as well as the augmentation radii were au- 
tomatically determined using the sphere-geometry opti- 
mization algorithm [l8[. Self-consistency was achieved 
by a highly efficient algorithm for convergence acceler- 
ation [191 ]. The Brillouin zone integrations were per- 
formed using the linear tetrahedron method with an in- 
creasing number of k-points in order to check the conver- 
gence with respect to the granularity of the k-point grid 
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The V2O3 surface was simulated by means of a fivefold 
superstructure of the canonical hexagonal unit cell. The 
latter comprises six formula units, i.e. six V layers along 
the (hexagonal) c axis. These layers are separated by O 
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FIG. 1: (Color online) Slab used for modeling the V2O3 sur- 
face. Large and small spheres represent the V and O atoms, 
respectively. 



layers, giving rise to VC>6 octahedra. For further details 
see, e.g., Refs. [U [22[ and the references therein. The 
lattice constants and fractional coordinates of the atoms 
were taken from Ref. [23|. The surface was then gener- 
ated by removing all atoms from one of the five hexagonal 
unit cells, see Fig. [TJ such that one of the resulting sur- 
faces terminates with an O layer. This vacuum interface 
is the subject of the subsequent discussion. 

In order to separate the charge redistribution due to 
the broken crystal symmetry and the changes in chemi- 
cal bonding at the vacuum interface from the structural 
relaxation effects at the V2O3 surface we will not include 
the latter in our considerations. In fact, structural re- 
laxation is restricted to a very narrow range of about 
five atomic layers as has been pointed out in a detailed 
study by Kresse et al. [HI, which also provides an elab- 
orate overview over the relevant literature. In contrast, 
we shall focus on the question whether the charge redis- 
tribution coming with the removal of bonding partners 
could affect the behavior of a much thicker slab of sev- 
eral 100-1000 A. Consequently, accounting for the near- 
surface structural modification is, on the one side, dis- 
pensable and, on the other side, disadvantageous for a 
systematic analysis of the decay of charge redistribution 
effects. 

Of course, in view of the strong electronic correlations 
present in V2O3, calculations beyond the standard LDA 
seem necessary. Indeed, while LDA fails to capture the 
metal-insulator transition of the bulk material, the lat- 
ter could be successfully described using LDA+DMFT 
[25l l26j . However, for the large supcrcell used in the 
present simulation of the surface the rather complex 



LDA+DMFT calculations would be beyond present day 
computational facilities. Even more important, the pre- 
vious LDA+DMFT studies revealed that the strong elec- 
tronic correlations come into play via the narrowing of 
the bands as calculated within the LDA, which occurs 
on going from the crystal structure of the metal to that 
of the insulator. Hence, there is a clear signature already 
in the LDA bands indicating the metal-insulator transi- 
tion. As a consequence, although LDA has turned out to 
be insufficient in describing the metal-insulator transition 
itself, it clearly signals its occurrence. For that reason, it 
is well justified to omit a full LDA+DMFT treatment as 
long as the LDA results are carefully interpreted. 

III. RESULTS 

Due to the octahedral O-coordination of the V atoms 
in V2O3, the V 3c? electronic bands (located around the 
Fermi level) are subject to a splitting in two subgroups. 
While the partially occupied t2 S manifold is found in the 
energy range from — 1.2eV to 1.4 eV, the unoccupied 
states extend from 1.7 eV to 3.7 eV. These two groups 
are clearly observed in Fig. [2l which depicts the site- 
projected V 3d density of states (DOS) as obtained for 
various atoms in the thin film. Atoms are numbered ac- 
cording to their distance from the surface with atom 1 
being located directly at the surface. We find that the 
shapes of the curves in Fig. [5] are closely related to each 
other. Only for atom 1, which is subject to the strongest 
interface effects, deviations are visible since half of the 
bonding partners are missing and the V-0 therefore is re- 
duced. In fact, the DOS curves of atoms 7 and 11, which 
are separated from the surface by about 18 A and 28 A, 
respectively, virtually resemble the bulk DOS shape. For 
this reason, the decay of the interface-induced relaxation 
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FIG. 2: (Color online) Site-projected partial V 3d DOS for 
a series of V sites in a V2O3 thin film. Numbering of the V 
atoms corresponds to their distance from the surface. 
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FIG. 3: (Color online) Partial V 3d and O 2p DOS (left) as well as V 3d symmetry components (right) for a VOe-octahedron 
in bulk V2O3 (top row), separated by ~ 28 A from the surface of a V2O3 thin film (center row), and at the surface of a V2O3 
thin film (bottom row). All DOS curves are per formula unit. 
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FIG. 4: Electronic band structure of bulk V2O3 (top left), for high-symmetry lines in the first Brillouin zone of the hexagonal 
lattice (top right), compared to a V2O3 thin film with bulk lattice constants (bottom left) and with an 1% elongated c lattice 
constant (bottom right). 



of the electronic states is found to be very fast. 

In order to establish further insight into the above- 
mentioned decay, we display in Fig. [3] V 3d and O 2p 
partial DOSs in a wider energy range in the left column as 
well as orbital resolved V 3d partial DOSs in the right col- 
umn. The nearly perfect splitting of the t2 S manifold into 
a\g and sub-groups traces back to the trigonal sym- 
metry of the corundum lattice (2f|. The partial DOS of 
V2O3 bulk (top row of Fig. [3]) is almost identical to that 
found for a VOg octahedron in the thin film some 28 A 
beneath the surface (center row of Fig. [3]). This holds 
for both the O 2p and the V 3d states. In contrast, the 
surface V06 octahedron clearly shows distinct electronic 
relaxation effects of the local electronic states (bottom 
row of Fig. [3]), which can be attributed to the modified 
chemical d-p bonding at the surface. In particular, for 
those V and O atoms, which miss a bonding partner due 
to the surface, we observe a reduction of the 3d and 2p 



occupation, respectively, of approximately 0.08 electrons. 
We point out again that in obtaining these results we did 
not take into account the V2O3 surface relaxation in or- 
der to distinguish between the electronic changes and the 
short-range structural distortions. Our results entail that 
the effect of the vacuum interface on the local electronic 
states, as reflected by the site-projected DOS, vanishes 
at about 15 A below the surface. Charge redistribution 
induced by the surface is therefore restricted to the same 
narrow surface layer, in which the structural relaxation 
is observed [24| ■ 

In a next step we turn to the discussion of k-resolved 
quantities and display the electronic bands in Fig. 2] The 
top row of Fig. 2] depicts the band structure of bulk V2O3 
as well as the non-primitive hexagonal Brillouin zone of 
the corundum lattice of V2O3, which here and in the fol- 
lowing is used for representing findings of band structure 
calculations. The bulk results serve as a reference for 
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comparison with the surface results given in the bottom 
left panel of Fig. [H 

Of course, the surface band structure comprises a lot 
more bands due to the simulation of the surface by the 
fivefold superccll of the surface slab as mentioned above. 
In this supercell, the increased number of different bands 
reflects the larger number of incquivalent atomic sites. 
V (and O) atoms with different distances to the vacuum 
are no longer equivalent from the crystallographic point 
of view. Furthermore, the dispersion along the high sym- 
metry line r— A is strongly suppressed due to the fact 
that this line corresponds to the direction perpendicular 
to the surface. Note that the construction of the sur- 
face slab in a fivefold supercell of the hexagonal V2O3 
cell implies that the length of the line r-A is reduced 
to a fifth and that the electronic bands are folded back. 
In order to enable a straightforward comparison with the 
bulk data, the line T-A is thus artificially stretched to the 
original length in the slab band structure. Having this 
fact in mind, we eventually find that the surface slab re- 
produces the characteristics of the bulk V2O3 bands. In 
particular, there is no sign of an additional state, not 
present in the bulk, which could be responsible for the 
anomalous behavior of thin films. 

It has been found for V2O3 thin films, that the a and c 
constants of the corundum lattice deviate slightly, i.e. less 
than 1%, from their bulk values 0,11]. In addition, we 
have shown that the influence of the a parameter change 
on the electronic structure is negligible 0. Here, we in- 
vestigate whether modifications of the c parameter may 
play a more important role by performing a surface slab 
calculation with the c constant increased by 1% as com- 
pared to the bulk. The results are depicted in the bottom 
right panel of Fig. dj which uses the same representation 
as before. Comparison with the bottom left panel of Fig. 
[4] shows that an elongation of the c axis does not lead 
to a qualitative change of the electronic states. The two 
band structures are, in fact, almost identical. Hence, we 
obtain that the experimentally observed alterations of 
the V2O3 lattice parameters in a thin film influence the 
electronic states only marginally, as long as they are not 
accompanied by a relaxation of the atomic positions. In- 
deed, a substantial lattice relaxation has been advocated 
to interpret transport measurements [2^ ]. 

To ensure that the thin film electronic structure fully 
resembles the bulk a final step is necessary. We still have 
to exclude differences in the contribution of the atomic 
orbitals to the different electronic states, since they would 
point at a modified chemical bonding. Orbitally weighted 
band structures are useful tools to address this question. 
Fig. [5] gives an example of a prototypical V atom, some 
28 A below the surface. Very similar results are obtained 
for the neighboring V sites. In Fig. El the bars added to 
each state have a length proportional to the contribution 
of the 3d t2 g orbitals of the atom under consideration. 
States with minor 3d t2 g admixture are not depicted for 
clarity. For this reason, some of the bands shown in the 
bottom left panel of Fig. 2] are missing in Fig. \5\ We 
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FIG. 5: Weighted electronic bands for a prototypical V atom 
in a V2O3 thin film, located m 28 A beneath the surface. The 
length of the bars shown for each state is proportional to the 
contributions of this atom. States with a contribution below 
a certain threshold are omitted for clarity. 

observe that the selected V site contributes to the 3d t 2g 
bands in the entire energy range from — 1.2eV to 1.4 eV 
and thus perfectly reflects the behavior of a bulk atom, 
which is true for any atom located at least ~ 15 A be- 
low the surface. Charge redistribution as induced by the 
vacuum interface therefore has to be excluded as a pos- 
sible source for the novel electronic properties of V2O3 
thin films. We finally mention that in the case of a het- 
erosystems, like LaA10 3 /SrTi0 3 [13, HI], different charge 
redistributions have been obtained in LDA and LDA+U 
calculations. However, for the present system we do not 
find such a dependence. 

In conclusion, our findings entail that neither the ob- 
served changes of the a and c lattice parameters in V2O3 
thin films nor the response of the electronic system to the 
presence of the vacuum interface can explain significant 
deviations from the electronic properties of bulk V2O3. 
Only next to the surface, i.e. in a narrow layer of at most 
15 A thickness, a charge redistribution is observed. The 
same behavior has been found in calculations for other 
transition metal surfaces [H, Hl[ . 

Below the 15 A surface layer the electronic states are 
hardly modified by the vacuum interface. LDA+DMFT 
calculations have shown that the narrowing of the V 3d 
t2 g group of bands seen in the plain LDA DOS (which 
traces back to structural alterations [32|) is the essen- 
tial criterion for obtaining a metal-insulator transition 
in bulk V2O3 [Hjlll]. Because such a narrowing is not 
observed in our results, they contradict a transition for 
V2O3 thin films of 100-1000 A thickness, as long as a re- 
striction of the structure relaxation to the surface layer is 
assumed. Reversely, only further changes of the atomic 
structure, probably in the entire thin film, can explain 
the experimental situation. A detailed study of the crys- 
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tal structure over the entire thin film is thus required to 
settle this issue. 
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